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TECHNICAL NOTE 3329 

SHOCKS IN HELICAL FLOWS THROUGH ANNULAR CASCADES OF STATOR BLADES 
By Robert Vasserman and Arthur V. Goldstein 


SUMMARY 

A method Is presented for calculating supersonic potential flows in 
annular cascades of blades by the method of characteristics. It Is found 
that helical flows may be adjoined by helical shocks of uniform strength; 
these constitute a considerable addition to the class of simple flows 
available for designing cascades of lifting blades. It was also found 
that by selection of the proper variables, the derivatives of the veloc- 
ity components, which occur in the characteristic equations, could be 
conibined into an exact differential. This form of the equation facili- 
tated computations. A flow and several cascade designs were computed. 


INTRODUCTION 

In studying the flow of air through a compressor, an analytic de- 
scription of the behavior of a particle on an "average" streamline or 
stream surface going through the congpressor Is quite useful. This type 
of analysis requires assumptions such as axial symmetry or flow on sur- 
faces of revolution. While they frequently serve well to get a good 
over-all picture of the flow, in certain regions of the compressor -these 
assumptions are generally untenable. These could be regions In which 
the major part of the loss Is occurring, so that such two-dimensional 
pictures would be inadequate either for design control or investigation 
of the cause of these losses. 

In order to construct compressors of maximum effectiveness and ef- 
ficiency, it is desirable to have a description of the internal velocity 
field. As a step toward this ultimate objective, effort has been applied 
to the less general problem of obtaining potential flows of a perfect 
nonviscous fluid satisfying boundary conditions required in compressors. 
Without assumptions of certain symmetries in the flow, such an analysis 
involves the solution of differential equations in three independent 
variables. Although several methods for obtaining general three- 
dimensional potential flows in compressors have been described (refs. 1 
and 2), the amount of numerical work required to get a specific solution, 
is formidable. 
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Thus, since on the one hand present two-dimensional methods axe 
frequently inadequate and on the other hand the existing three- 
dimensional methods are impractical with standard computing techniques, 
an attempt was made to find simple classes of flows, or close approxi- 
mations thereto, which could he used with the houndary conditions re- 
quired in a compressor. In particular, interest was focused on obtain- 
ing supersonic flows for the inlet region of a diffusing cascade of 
stator blades. From the class of uniform flows, a flow can he construct- 
ed through a straight cascade of lifting blades by joining such flows 
with plane shock surfaces which originate on the blade surfaces. The 
velocity potential is, in this case, a linear function of the cartesian 
coordinates. For an annular cascade, a free-vortex upstream flow results 
in shock surfaces which must be oriented for uniform stock strength in 
order for a velocity potential to exist . The simplest class of such sur- 
faces are those linear in 9 and z, which are the polar angle and axial 
distance, respectively, of the cylindrical coordinate system. These 
shock surfaces impose certain boundary conditions on the velocity com- 
ponents as functions of r (distance from polar axis) for the flow on the 
downstream side. It was found that the solutions which satisfy these 
boundary conditions are helical flows which have velocity components ex- 
pressible as functions of two independent variables. The helical flows 
for annular cascades are therefore the counterparts of the uniform flows 
for straight cascades in that they represent the simplest class which con- 
tains the free-vortex upstream flow and which may be adjoined by uniform 
shocks, while the uniform (constant) flows represent the simplest clasB 
which contains the constant upstream flow and which may be joined by 
shocks . 

Because the upstream flow and the shock-surface normals involve 
only the variable r, a velocity potential of the form <p =bz + c0 + g(r) 
(where b and c are constants) might be tried; that is, the simplest 
class of potential flows with a radial component. However, the boundary 
conditions imposed by a shock of uniform strength give functional rela- 
tions between the velocity components that cannot, in general, be satis- 
fied by a flow of this class. Another approach was to try to adequately 
approximate the potential function in terms of arbitrary functions of 
r and then satisfy the boundary conditions exactly. This procedure re- 
sulted in inadequate approximations to the potential -flow equation. In 
the end, no potential flow or approximate potential flow, expressible in 
terms of functions of single variables, could be found which satisfied 
the required boundary conditions or approximations thereof. 

Thus, while the three-dimensional problem is reduced without further 
assumption to a two-dimensional one, further reduction seems unlikely. 
Consequently, the problem was solved by the method of characteristics in 
two dimensions. In the mathematical sense, the flow is two-dimensional 
in that the flow problem involves two Independent coordinates. In the 
physical sense, the flew Is three-dimensional in that the velocity vector 
fields are not identical for a set of surfaces as in the cartesian flow. 
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SYMBOLS 

The following symbols are used in this report: (All velocities are 

referred to the stagnation sonic speed as the unit of velocity.) 

A,B constant coefficients in shock equation P = 0 and in coor- 

dinate transformation 

a local velocity of sound 

C,D constant coefficients in coordinate transformation 

F,G functions of coordinates which define shock surfaces 

f function of r in shock equation F = 0 

h component of stream function w, h = h(r,s) = w + t 


Jacobian of transformation of coordinates, J = f? 

<Hr,s,t) 


1 

sin a 


cos p 


sin 2 a + 


sin 0 ^5— 

v r l 78 1 r(Vs) 2 


Zf\ 


rWp , » [ a 2 sin p 2T t^l , (r-DvjW* 'l 

_r 2 (Vs) 2 Ur 2 (Vs) 3 2T 2 r|a 2 C ° S J 

M Mach vector in Vr,Vs-plane (M = u/a) 

M n Mach number normal to shock surface 

R r ,R 0 ,Rk base vectors along helical coor dinat e curves, 

F _ VsxTt , 

** ~ 7r.VsxVt' etc 

r cylindrical and helical coordinate (distance normal to axis 

of coordinate system) 

s helical coordinate, A0 + Bz 

T corrected stagnation sonic speed, T =A/1 - v 2 /R 2 

7 2 t^ t 


t 


helical coordinate, C 9 + Dz 
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U 

U 

u 

V 

V 


magnitude of U 

velocity component in Vr, Vs -plane 
integral of streamlines independent of t 
magnitude of V 

velocity vector, V = V r Vr + VgrVS + V z Vz 

physical velocity components in cylindrical coordinate system 
v t’ v B> v t V-R,,, V-Rt, respectively 




V r , V s , v 6 

W 

W 


z 


V*Vr, V*Vs, V*Vt, respectively 

ratio of U to corrected stagnation sonic speed, \T = u/T 
integral of streamlines linear in t 

axial distance of cylindrical coordinates (distance parallel 
to axis of coordinate system) 


a Mach, angle, a = sin" 1 ^ = sin" 1 i = sin" 1 

u M 

3 flow angle in Vr, Vs -plane measured from 

Vs, (3 = tan" 1 / - — “ 7 ~\ 

\| Vs jv 37 / 

T ratio of specific heats 

A AD - BC 

S -2 1 - M| 

r+1 M n 

e angle of deflection of flow in Vr,Vs-plane through shock 

wave, e = p 2 - or e s p 4 - p 3 

£ + ,£ slopes of characteristic curves in r, s-surface 

0 polar angle of cylindrical coordinates 


a 

TW 

Vr toward 
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/ cEW 
cot a 


g,T] characteristic coordinates 

p density of gas (relatiye to stagnation value) 

x leading-edge vector 

$ potential function 

Subscripts: 

refer to fig. 1 

d,e,f 


h hub 

0 leading edge of blade 

1 upstream side of F = 0 

2 downstream side of F = 0 

3 upstream side of G = 0 

4 downstream side of G = 0 


CONSTRUCTION OF SUPERSONIC DIFFUSOR FLOWS 
Passage of Flow Through First Shock Surface 

The flow upstream of the cascade is assumed to be a supersonic free- 
vortex flow; that is, V r = 0, rVg = constant, V z = constant, and the hub 

and the casing are circular cylinders. Also, the suction surface of the 
blade Is initially a stream surface of the upstream flow. When the 
leading-edge wedge angle is not zero, an attached shock comes off the 
leading edge on the pressure side of the blade. This shock will Inter- 
sect the hub, casing, and suction surface of the adjacent blade In cer- 
tain curves. 

In order that the flow downstream of the shock be a potential flow, 
the change of entropy across the shock must be constant and therefore 
the Mach number normal to the shock must be constant. If the shock sur- 
face is given by F(r,0,z) = 0, then F must satisfy 
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Vn op 

^ * ]Wf = Mn,l = constant (l) 

on F = 0. Because the upstream velocity components are functions of 
r alone, the coefficients of the differential equation (1) for F do 
not involve the coordinates 9 and z. A function F of the form 

F = AS + Bz + f (r) 

vill therefore satisfy equation (l), which reduces to 



where s is the quantity AS + Bz, and a is the sonic speed given By 

a 2 = 1 - V 2 (when the stagnation sonic speed is taken to he the unit 
of velocity). Thus, the class of shock surfaces obtained is 

f(r) + s=0 (3) 

in which f is the solution of equation (2) . If it is required that 
the casing he a cylindrical surface of revolution, a reflected shock sur- 
face G = 0 results which intersects the surface F = 0 at the casing 
and which is so oriented that no radial deflection to the flow will re- 
sult from passage through both shocks. The shock surface G = 0 will 
intersect the huh and the blade surfaces downstream of F = 0. The 
shapes of the huh and suction surfaces downstream of F = 0 are ob- 
tained by the condition that the shock F = 0 is not reflected at these 
surfaces . 


CVJ 


a 


to 


Flow Field Downstream of First Shock 

After calculation of the flow through the first shock (F = 0), the 
downstream flow field must be constructed.. Helical flows will be shown 
to satisfy the boundary conditions on the shock, thus reducing the prob- 
lem from a three-dimensional to a two-dimensional flow not restricted 
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"by the usual symmetries required of two-dimensional flows in cartesian, 
cylindrical, or spherical coordinates. The problem will he further re- 
duced to ordinary differential equations by the introduction of charac- 
teristic coordinates, and then further simplified by the introduction 
of the Erandtl-Meyer angle variable. 


To obtain the flow downstream of F = 0, consider the boundary con- 
ditions imposed by this surface. For the velocity downstream of F = 0, 

? 2 = Vl - (4) 

where 


8 -2 1 - 
1 r+1 


and the signs of df/dr and the coefficients A and B are chosen to 
give VF*V 1 > 0 and V r ^ 2 < 0. In terms of a helical coordinate sys- 


tem given by 


r = r 

s = A9 + Bz ’• 
t = C0 + Dz 


(5) 


the boundary conditions by equation (4) are 


_ _ , a-i &-i — a-i &-i f * 

V*r = ^ f ‘ = -[VT]- 




— — 8-1 ~\ a-i &-i 

y 2 .E B = V r E B - py = S (DrT 9 - OT Z ) - 


( 6 ) 


Vg*R^. = ^ (- BrVg + AV Z ) = constant 

where A s AD - BC, = (VsxVt)/(Vr*VsxVfc), and so forth. 


If there is a potential flow 
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with components 


v r = Rp-V 


= Rp-VS = 


S$ 

5? 


V S - V* = Sy 



and if it is of the degenerate type having Sv_j/St = 0 (helical flow) 
where i = r,s,t, then 

Sv t Sv s n 

5T" ** 5t~ = 0 
Sv + Sv_ 

and v^. is therefore a constant. The houndary conditions (6) fulfill 

this requirement. Consequently, the equations for the flow downstream 
of F = 0 are written in r,s,t coordinates, and solutions are sought 
for which Sv^/St = 0. 

The continuity equation is 

- V- (pY) = V-Y + V--^. _ o 
P ' P 


Because 


then 


p = constant 



1 


r-i 

2 



V-Y - Y-VV 2 = 0 
2a^ 

The formula for the divergence (ref. 3), in the case S/St = 0, reduces 
to 
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W 

w 

CO 


7 - v - 7 [5 ( J7r - V > + h (J7s-,f) ] 


where J is the Jacobian of the transformation : 

T = _ 1 _ r 

d(r,s,t) 7r*7sx7t A 

The second term of the continuity equation is 


(7) 


- ^2 5 (V* + ^ + V‘) + ^ ^ (v^ + ^ + T t T*)_ 


CO . 

I 

H 

O 


where 


v 1- = 7r*V = 7r* (v r 7r + v g 7s + v^.7fc) = v T 
% 

V s s Vs*V = v s (Vs) 2 + v.j.(7s«7t) 
v* = 7t*V = v g (7s»7t) + v 1 .(7t') 2 




( 8 ) 


J 


Because v.j. is constant, the equation of continuity may be reduced to 


[■ 


(v*T - a 2 


r~- + Z'Fv 5 + (v 6 ) 2 - a 2 (7s) 2 - — a 2 +r 2 (V*70) 2 


= 0 

(9) 


where derivatives of v 37 and V s have been eliminated, and the equa- 
tion for zero vorticity 


Sv r dv s 
3T - " 5r~ 


= 0 


( 10 ) 


has also been employed. 

If characteristic coordinates £ and tj are introduced, then the 
equations for r and s in terms of £ and tj are 


9s y dr 9s * dr 

3£ = 5if 


(ii) 
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where £ satisfies the quadratic equation 

[a 2 - (v 1 ) 2 ] i 2 + 2v r v s £ + [a 2 (Vs) 2 - (v 8 ) 2 ] = 0 
(ref. 4) with two roots 


__ -v r y s ^ a}?sl (v 2 *) 2 + ( V s /! Vs 1 ) 2 - a 2 


, 2 - ( v ^) 2 


( 12 ) 


(13) 


The characteristics flow equations are simplified somewhat "by uti- 
lization of the following relation from equation (12 ) : 




a 2 (?s) 2 - (v s ) 2 
a 2 - 


and the form which results is 


dv r 

w 


+ 


dv 

(O w 


v 2 * 

r 


a 2 + r 2 (Dv 8 - Bv^/A 2 
a 2 - (^y 2 


hr 




0 


> 


dv r 


+ 


, Sv s . v 2- Fa 2 + r 2 (Dv 8 - Bv t ) 2 /A 2 "|&r 

t 5 + J 3 T r [ a 2 -(^>2 


0 


J 


(14) 


Equations (ll) and (14) are to he solved in a region downstream of 
F = 0, with the boundary conditions on F = 0 given by equations (6) . 
Present experience by the authors in the use of these equations indi- 
cates that, when the trapezoidal rule of integration is used with finite 
differences, an integral converges more rapidly when the velocity incre- 
ments given by the first two terms are combined into an exact differen- 
tial. This can be accomplished by introducing the Prandtl-Meyer expan- 
sion angle v and the direction angle 0 as the variables defining the 
flow. The velocity vector V may be expressed in terms of the set of 
unit orthogonal vectors Vr, Vs/|Vs|, and Rj./| Rj. J as 


v 8 Vs y t ^t 

RFI l^TT^T 


V = T^vr + 
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oi 

CO 

CO 


The velocity components v 1 * and V s /] Vsj give a vector U, where 


U = u|c 


cos 0 Vr + sin 0 



Therefore, 

v 37 = U cos 0, v s /|vs| = U sin 0 

Since v r = v 1 " and 

V s = Vs-V = Vs*(v r Vr + v s Vs + v^Vt) = v s (Vs) 2 + vj-(Vs-Vt) 


then 


v r = U cos 0 


_ U sin 0 Vs.yt 

Vs | vs | " H (vb ) 2 

The characteristic slopes ^^/jvsj involve not only U and 0, 

hut also a 2 = 1 - (U 2 + v 2 /r|) and therefore the coordinate r hy 

means of R 2 . However, if a modified stagnation sonic speed T is in- 
troduced where T 2 = a 2 + U 2 and the velocity is referred to this 

speed, then an exact differential may he constructed. Let W he the 
ratio of velocity U to the modified stagnation sonic speed 


When the Mach angle a is introduced, where 
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there are obtained 


5± , _ , 

jfTj- = tan CP + a ) 


and 


ds = | Vs j tan (p + a) dr 

in place of equations (13) and (ll) . With the definition 

v mf cot 

the flow equations (14 ) become 


(13a) 


(11a) 


where 


dv ± dp 



(14a) 


% sin a ) cos 


cos p [* . o , / A . . AB v t\ 2 

(p + a) L \t Vs r(Vs) 2 u / 


sin (P ±. a) 


a 2 2v t AM 1 (r-l)v 2 B 2 r 2 

sin p — — ^ + o , „ — cos a ' 


r 2 (Vs) 


Ur 2 (vs) 3 J 2T 2 B^A^ 


For the method of characteristics to he useable, a must he real, and 
therefore U > a. That is, the velocity component in the Vr, Vs-surface 
must he supersonic regardless of the magnitude of v^.. 


For trapezoidal integration, the following formulas are used to 
calculate conditions at point c (fig. 1(a)), where c is connected with 
point a hy a curve of tj = constant and with point h hy a curve 
£ = constant: 

v c “ v a Pc “ Pa £^+^c,a (l°S r c “ = ® ''j 


s c - s a = [l? s i tan CP - a )] c ,a C r c - r a) 
v c ~ v b " Pc + Pb + t K J C ,b (log r c - log 1^=0 

s c - ®b = [i Vs J tan (P + a)] c ^ (r c - i^) 


( 15 ) 
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where the square "brackets indicate the average of values designated ■by- 
subscripts . The four equations will determine the -four quantities v c , 

0 C , s c , r c . Iteration is required for accurate values of the Z*s. 


Continuation of Flow Field Through Reflected Shock 

To continue the solution downstream, the condition that the casing 
he a cylindrical surface (v r = 0) is applied. The calculations up to 

this point cover the region ABC (fig. 1(h))-, which is hounded hy the 
first shock and the two characteristic curves £ = constant (AB) and 
T] = constant (BC). The required reflected shock is again assumed to 
satisfy equation (l); that is, if the shock is given hy G(r,s,t) = 0, 
then G satisfies 


M • 


VG 
1 VGl 


= ^,3 = 


constant 


(16) 


on the surface G = 0. However, since the coordinate t is not explic- 
itly involved in any of the coefficients of the derivatives of G in 
the homogeneous differential equation (16), it is possible to find a 
function G(r,s) of r and s only, which will satisfy the condition 
(16). Furthermore, since on the boundary 


that is, 

v r, 4 = ®r’^4 

v s, 4 - VV 4 


v 4 ~ V 3 “ a 3 B 3 


v r,3 " a 3 8 3 H/M = 


a function of r or s 


* dG 
v s,3 - a 3 6 3 SF 


/|7G| = 


a function of r or s 


v t,4 = V v 4 = % 3 = constant 


the boundary conditions are compatible with the assumption of a helical 
flow in the field downstream of the shock surface G = 0. That is, v r 

and v g are functions of - r and s alone, and v^. is constant. As 

before, the problem is a two-dimensional one in the r, s-plane (contain- 
ing Vr, Vs). 
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A typical situation which arises in the calculation is shown sche- 
matically in figure 1(h) . The velocities are supposed known on a curve 
def which crosses the shock wave; the problem is to extend the solution 
first to the point P and subsequently to a new region containing the 
point P and the extended portion eP of the shock curve. The normal 
Mach number may be regarded as known, since it is constant and equal to 
the value determined at the tip (point C, fig. 1(b)) from the velocity 
in the Vr, Vs-plane and the deflections of the flow at the shock. 


In addition to the shock equations, there are also available equa- 
tions relating the velocity parameter v, the direction angle p, and 
the position of the point P with those of the points d and f, which 
are located on the characteristics curve tj = constant. For trapezoidal 
integration these equations are 


v 3 " v d " ^3 + ^d + ^ K -^3,d( log r 3 " log r d) = 0 

/ 

v 4 - v f " P 4 + Pf + £ K J 4 f ( log r 4 - lo S r f ) = 0 


(17) 


s 3 ' s d = [l Vs I tan (p + a)J^ d (r 3 - r d ) 

S 4 “ S f = [ I 7s 1 tan + a ^4,f ( r 4 " r f ^ 


(18) 


where the subscripts 3 and 4 refer to the point P upstream and 
downstream of the shock surface, respectively. The double subscript and 
square brackets refer to an average between the points. It is necessary 
to locate d and f in order to obtain by interpolation the values of 
the quantities required. By utilizing the slope of the shock and the 
curve tj = constant, there are obtained the intersections 


log r d = log r g 


[r|vs|tan(p+a)] 3 d -[r|Vsltan(p+0w)] 3 ^ e 
+ [r jVs |tan(p+a)j ^ [r |V s | tan(P-a)] d ^ e ^ og r 3 " og 


*e> 


*\ 




log r_= log 


[r| ysj tan(p+a)] 4jf -[r j 7s [ tan(p+9w) j 5j e 
T S + [r j Vs|tan(p+cc.) ] 4 ^-[r ) Vs Jtan(P-a)] g ^ f ° g 


r 3 ” log r e ) 
(19) 


where 0w is the wave angle of the shock and tan (p + 0w) is the slope 
of the shock wave. By subtraction of the characteristic equations (17 ) 
and utilization of 


P 4 = P 3 + e 


( 20 ) 
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where e is the deflection angle of the flow passing through the shock, 
there results 

v 3 - v 4 + s = ( v d - 0 a ) - (v f - 0 f ) + tKj 3 ^(log r d - log r 3 ) - 

C^_3^f ( lo S r f - log r 4 ) (21) 

Equation (21) may he solved in conjunction with the condition of 
constant normal Mach numbers and the shock relations hy the following 
process: (l) Assume values of Mgj then find values for v 3 and ctg 

from supersonic-flow tables. (2) Calculate 0w = sin -1 ^ 3/M3). (3) 

From shock tables and values of M 3 and 0w, find values for e, M 4 , 
and v 4 . (4) The correct solution is that for which Vg - V4 + e is 

equal to the value calculated by equation (21) . 

By iteration a convergence might be reached in which equation (21) 
and the shock relations are first solved only roughly for values of Vg, 

V4, 0 3 , and P4. Solutions are then obtained for r d and r f , followed 

by interpolation for the required variables at those points. Thp proc- 
ess is then repeated with refined values of the coefficients. However, 
this process will often fail to converge for weak shocks because the 
condition of constant normal Mach number and the condition on 
v 3 - v 4 + 6 are nearly the same. Therefore, the solution is nearly 

indeter m i n ate and will fluctuate widely from one approximation to the 
next. That these two conditions are not suitable for determination of 
the shock wave may be seen from the following considerations. When the 
shock is weak, Vg - V4 « c, and equation (21) gives a value for 

v 3 - v 4 + 6 « 2s. Also, a condition of constant normal Mach number 

corresponds to constant pressure ratio. Graphs of shock-wave so lu tions 
will show that, for assigned values of the pressure ratio and a, the 
incoming Mach number is nearly indeterminate. 

Since for weak shocks the solution is practically indeterminate, an 
alternate procedure was used, which involved the assumption of a value 
for one parameter such as P 4> The value of V4 may then be determined 

from the second equation (17 ). By using the approximation V3 = V4 + e 

in equation (21), a value is then determined for &} this value is sub- 
stituted into equation (20) to determine pg and v 3 = \ + *• Tables 

are used to find Mg from v 3 , and the wave angle is then calculated 
from 0w = sin“^-(Mg/M I1 g) . The results should be checked for consist- 
ency with the shock tables. If the agreement is not close enough, a new 
value for P4 is assumed. 
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It has been found expedient in the solution to fix a value for 
and to permit S3 and the points d and f to fluctuate from one 
iteration to the next, because the locations of d and f vary only 
slightly thus permitting rapid convergence and because the size of the 
network is controlled. 

When the solution for the shock is started at the casing, the situ- 
ation is somewhat different in that p = 90°, but a value for v is not 
available. However, if the point P is chosen close enough to C that 
variation in parameters may be neglected from C to the point on the 
casing corresponding to f, no great increase in computation will result 
because the size of the interval on the shock curve can be rapidly in- 
creased after the solution is begun. Such a procedure is valid, because 
it can be shown that it results in a solution continuous at the point C. 


BLADE AND HUB SUBFACES 

After the flow field has been obtained, the hub and blade surfaces 
may be computed by using the property that they are stream surfaces and 
therefore contain the velocity vector V. Because the equation of con- 
tinuity 

V*(pV) * 0 

is satisfied, there exist two independent stream functions u,w which 
satisfy the equations 

Vu-V = Vw-V = 0 

These functions may be "identified by exa m ini n g the expanded form of the 
continuity equation 

|r (pjv*) + ^ CpJv*) + (pJv*) = o 

where J (= r/A = (Vr-VsxVt) -1 ) is the Jacobian of the transformation of 
coordinates. Because all the quantities pJvi are functions of r and 
s only, then the last derivative is zero and the remainder of the equa- 
tion indicates the existence of a potential function u(r,s), which is 
the stream function of the flow . Then, 


or 


dn = pj^dr - v’mIs) 


or 
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Vu = pJ(v®Vr - v^Vs) = pV*R{. 

The surface u = constant is therefore the solution of 


dr ds 

" 7 ® 


With the first solution assumed known, then the second stream function 
w can "be found such that 


pV = VuxVw 

5 For an arbitrary path element dR, 

p dRxV = Vu dR*Vw - Vw dR*Vu = Vu dw - Vw du 

whereas for a path element on the surface u =. constant, there results 

h dw = = PA‘<iR*V _ _ dR-A*Y 

° A*Vu pA*V*Rf. " ifAxv 

where A is an arbitrary vector such that A- Vu 0. If, for example, 
A = Rp, then 

v* v t 

dw = - dt + — ds = - dt + — dr 

V s v 1 

where, in the course of the integration, u iB to he regarded as con- 
stant. In integrated form. 


w + t = h(u,s) = — ds 

J u v 8 

The coefficients of integration may he exchanged for those computed in 
the field, so that 

u = | J' p Ur ( | Vs | sin p dr - cos (5 ds) 

Since u may he calculated hy integration on any path, the integrations 
are made on characteristic curves; these curves were chosen because the 
data have heen computed on them. Therefore, 

, ds = jVs] tan (p ± a) dr 
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and 


n = ± 


;/ p 


Urjv s | 


cos 


sin a dr 
(P T a) 


-*x 


Y+l 

a"*" - '*' r|vs| 


cos 


TFT 


dr 

sr 


The upper sign is used for integration on curves of t) = constant and 
the lower, for curves of g = constant. The sonic speed is given by 



¥« 2 ) 


1 + M 2 


where 

T s Vl - ^ v|/R| 


and M = l/sin a. That is, a/T is taken as the ratio of sonic speed to 
stagnation sonic speed, which is consistent with the values of v, a, or 
M. 


After contours of u = constant are established, then at appro- 
priate points on the contours r, a, and a are determined and the 
following integration is made: 


h(u,s) 




ds/ Vs 
U sin p 



ds 


with 


U = a/sin a = Tff = — ========== 

Vl + ^M 2 

The leading edge of the blade is defined by a curve that lies in 
the surface F = 0 and on which r assumes a set of values r^<r <r^., 

where r^ is the inn er wall radius. Consequently, r may be used as a 
parameter for the curve 


s o = s o( r o) 
*0 “ t o( r o) 
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where the subscript 0 indicates the leading edge. However, the con- 
dition that F = 0 establishes a condition on Sq, and the leading 

edge is given by 


s 0 = - f ( r o> 

*0 = *0^0) 

When h(r, s) is calculated, the integration on surfaces of constant u 
may be started from the curve F = 0 with a boundary value of h = 0, 
so that h(r, - f(r)) = 0, and consequently hg = 0 and 

w 0 = h 0 - t 0 = - t 0 (r 0 ). A sequence of values of r Q will result in a 
sequence for w 0 , Sq and therefore Uq, so that the the following func- 
tional relation is given between w and u, which defines the stream 
surface in question: 


w = - *0^0) 

u = u(r 0 ,s 0 ) = Uo(r 0 , - f(r Q )) 

A section of this surface at constant t is obtained by assuming 
values of h, which in turn define w « h - t. The stream-surface rela- 
tion between w and u then gives the 'correspondence between h and 
u, thus determining the desired curve. A section at constant r is ob- 
tained by assuming values for s, which then define h, u, w, and 
t = h - w. A section at z = constant is determined by substituting 
the value of t in terms of h and w(u) 

t = h - w(u) 

in the following relation between z, s, and t: 

Cs - At = (BC - AD)z 


Then, 


Cs - Ah + Aw(u) = (BC - AD)z 
expresses a relation between r and s that, with 

A9 - s - Bz 


determines the section. 
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If A = 0 and z is constant, then s = constant, 
t = 0 = h(r,0) - w(u) gives the desired relation between r and 0 
for constant z . 

Finally, the factors involved in the selection of the function 
t 0 (r 0 ) are considered. One factor of importance in fabrication is the 

wedge angle of the blade as measured between the intersections of the 
blade surfaces in the region near the leading edge with plane elements 
normal to the leading edge. These angles are obtainable directly as the 
deflection angles for oblique shocks where the gas velocities are com- 
ponents normal to the leading edge. If the leading-edge vector is given 

"by 


* = Br + *s§ + Rt§ 

the condition that t*VF = 0 determines ds/dr *= - f'(r). Then, 

dt 


T = Ej. - f * (r)R s + R,. 


dr 


The magnitude of the velocity required is then |txV(/|t|, and the ef- 


fective Mach number and wave angle are 


ItxvI 

a I* | 


V |w| )’ 


and Bin" 

respectively. The deflection angle can then be evaluated from standard 
tables as a function dt/dr and from the initia l shock surface. 


EXAMPLE 

As an example, a flow field was calculated for the following condi- 
tions: The free-vortex upstream flow was determined by the assumption 

of a Mach number of 2.0 at an angle of 45° with respect to the axial 
direction and lying in a cylindrical surface r = 1.0. The curve rep- 
resenting the intersection of the initial shock with the surface 
r = 1.0 was determined by setting s = 0. The shock strength was also 
assumed by setting = 1.0927. 

With no loss in generality, C = 0, D = 1. The for mu la for K^. 
reduces to 


K ± = ± tan (p T a) 
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In constructing the network of points not on a shock, the first 

calculation gave accurate values of v , 0 , r , 0 when the assunrp- 

c c c c 

tions that [K + ] c ^ a * [K + ] a and £K_3 c ^- b * [K.]^ were used. However, 
the coefficients [K + ] c , [K_] c , tan (0 C - a c ), and so forth, were calcu- 
lated and one iteration executed as a check. This calculation was based 
on a net with an increment in r of approximately 0.05 "between adjacent 
computed points of a characteristic curve. "When a point for extension 
of the shock curve G = 0 was calculated (F = 0 was calculated "by a 
direct integration) , the first iteration gave an accurate answer and a 
second was required for a check. If good estimates are made on the 
variation of the coefficients, the first iteration could serve as a 
check on the first integration. 

The net of characteristic curves and the initial and reflected 
shocks are drawn in figure 2 in the surface normal to the diffuser axis 
and viewed in the direction of flow. The Mach number contours resulting 
from the velocity component normal to the surface downward through the 
sheet and the calculated r,0 components are also shown. These contours 
are, of course, discontinuous at the shocks. 

Stream functions are shown in figure 3. The values of u have "been 
divided "by u^, which is the value of u for the entire flow in the 

region 0.7 < r < 1.0 upstream of the first shock, in order to indicate 
more clearly the equal increments into which the flow has "been divided. 
Similarly, the contours of h = constant are shown with values of h 
divided "by a displacement of z equal to 0.06556 times the tip radius. 

If a "blade is assumed to pass through the shock surface F = 0 at 
z = 0, then the h contours are "blade pressure-surface sections at con- 
stant z, spaced at axial-distance increments of 0.06556 times the tip 
radius. If 20 blades are assumed, then at r = 1.0 the preceding blade 
is encountered at 0 = - 18°, which corresponds to h/o. 06556 = - 4.8. 
These sections were computed and then shifted 18° to indicate the pres- 
sure and suction surface of a blade at the same z -location rather than 
the channel between two blades. The hub surface to the right indicates 
the hub at the suction surface. 

A different shape for the blade can be obtained by assuming 
w = w(u) instead of w = 0. For example, if the pressure surface of the 
blade is to contain a radial line, then, when z = 0, 0 must be zero. 
Therefore, on the initial curve, w(u) is established by 

w(u) = h(u,0) - z = h(u,0) 

From this equation, h can be determined for each pair of values u,z, 
and consequently the value of 0 can be determined from the value of 
h(u,0) and of u. Curves of this type are shown in figure 4, where the 
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■blade sections are taken at equally spaced values of z. The suction 
side of the "blade is computed as for the previous example (fig. 3), with 
the assumption of 20 blades. 


EVALUATION OF METHOD 

Calculation of helical flows by characteristics requires more time 
than for cartesian flows, but the time involved is not prohibitive. The 
main advantage of the method is that it gives rigorously a more general 
class of flows than have been heretofore available by two-dimensional 
methods. It is now possible in a larger variety of circumstances to 
give a q ualit ative description of actual flows in three-dimensional ma- 
chines and to broaden the approach to the design problem. As a design 
method, it has the disadvantage of not permitting the description of the 
blade shape, hub, and casing surfaces in advance of the calculation. 
However, once a flow is computed, a variety of blade shapes may be com- 
puted, depending upon the assumed spacing and leadi n g-edge orientation. 
These blade shapes Include swept and tilted blades, thus permitting a 
new degree of freedom to the designer. 

A restriction on shock-wave orientation for cascade design is evi- 
dent on comparison with the class of uniform flows. TheBe uniform flows 
may be joined by plane shock surfaces of any orientation provided that 
the normal. Mach number is supersonic. The helical flows, on the other 
hand, when joined by shocks of uniform strength, generally give rise to 
flows with a variable helical covariant velocity component v s in addi- 
tion to the radial component. Shock surfaces are therefore limi ted to 
the class containing the vectors R^. normal to both Vs and Vr. That- 

is, the coordinate s must not be changed as long as v s is not con- 
stant. When v_ is reduced to a constant, then a new class of shocks 

may be used and the variable s may be modified. In the case of the 
nmniiar cascade this does not represent a limitation in design procedure, 
because the uniform flows are inapplicable and the helical flows are 
required. 


Lewis Flight Propulsion Laboratory 

National Advisory Co mmit tee for Aeronautics 
Cleveland, Ohio, September 14, 1954 
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Figure 2. - Characteristic curves, Shock curves, and constant Mach contours. 
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Figure 3. - Stream functions u, h and axial sections of blades 18° apart pass- 
ing through z = 0 on Initial shock. 
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Figure 4. - Axial sections of "blades 18° apart and passing through 8 = 0, z « 0. 
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